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Abstract 

In this paper we describe a wide class of non-Volterra quadratic stochastic operators using 
N. Ganikhadjaev's construction of quadratic stochastic operators. By the construction these 
operators depend on a probability measure being defined on the set of all configurations 
which are given on a graph G. We show that if [i is the product of probability measures being 
defined on each maximal connected subgraphs of G then corresponding non-Volterra operator 
can be reduced to m number (where m is the number of maximal connected subgraphs of G) 
of Volterra operators defined on the maximal connected subgraphs. Our result allows to study 
a wide class of non-Volterra operators in the framework of the well known theory of Volterra 
quadratic stochastic operators. 

Keywords: Quadratic stochastic operator, Volterra and Non-Volterra operators, simplex, 
product measure. 

1 Introduction 

It is well known that the principles of biological inheritance, initiated by Mendel allow an exact 
mathematical formulation. For this reason classical genetics can be regarded as a mathematical 
discipline. The fundamental investigations on these problems were carried out by S.N.Bernstein 
[2] and Reiersol [13]. 

The main mathematical apparatus of such investigations, to our knowledge, is the theory of 

quadratic stochastic operators. Such operator frequently arise in many models of mathematical 

genetics [1-11], [13], [15], [16]. 
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The quadratic stochastic operator (QSO) is a mapping of the simplex 

n 

S"- 1 = {x = (x 1 ,...,x n ) G R n : Xl > 0,J2^ = 1} 

i=i 

into itself, of the form 

n 

v : x 'k = Pij,kXiXj, {k = 1, ...,n), (1) 
where p^fc are coefficients of heredity and 

n 

Pij,k>0, ^Pij,k = l, (i,j,k = l,...,n). (2) 

k=l 

Note that each element x G S 71 " 1 is a probability distribution on E = {1, ...,n}. 

The population evolves by starting from an arbitrary state (probability distribution on E) 
x G S*™" 1 then passing to the state Vx (in the next "generation"), then to the state V 2 x, and 
so on. 

For a given x< ' G S""" 1 the trajectory {x^}, Z = 0, 1, 2, ... of QSO (1) is defined by x<- l+ V = 
V(x®), where I = 0, 1, 2, ... One of the main problems in mathematical biology consists in the 
study of the asymptotical behavior of the trajectories. This problem was fully solved [1-7] for 
Volterra QSO which is defined by (1), (2) and the additional assumption 

p ijik = 0, if k (£ (3) 

The biological treatment of condition (3) is rather clear: the offspring repeats the genotype of 
one of its parents. 

In paper [6] the general form of Volterra QSO V : x = (xi, ...,x n ) G S" 1 ^ 1 — > V(x) = x' = 
(x[, ...,x' n ) G S" 1 ^ 1 is given: 

x' k = xAl + Y^ a kiXi ) , (4) 
^ i=l ' 

where a k i = 1Vik,k — 1 for i 7^ fc and a k k = 0. Moreover a k i = —a-ik and \a k i\ ^ 1- 

In papers [6], [7] the theory of QSO (4) was developed using theory of the Lyapunov functions 
and tournaments. But non- Volterra QSOs (i.e. which do not satisfy the condition (3)) were 
not in completely studied. Because there is no any general theory which can be applied for 
investigation of non- Volterra operators. To the best of our knowledge, there are few papers 
devoted to such operators (see e.g. [1], [4]). 

In papers [3], [5] a constructive description of QSOs is given. This construction depends on 
a probability measure ji and cardinality of a set of cells (configurations) which can be finite or 
continual. In this paper we describe QSOs using the construction of QSO for the general finite 
graph and probability measure /j, which is product of measures defined on maximal subgraphs 
of the graph. We show that if /x is given by the product of the probability measures then corre- 
sponding non- Volterra operator can be reduced to m number (where m is the number of maximal 
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connected subgraphs) of Volterra operators defined on the maximal connected subgraphs. Then 
using the well known theory of Volterra operators we study the asymptotical behavior of the 
trajectories of the non- Volterra operators. 



2 Construction of QSO 

Each quadratic operator V can be uniquely defined by a cubic matrix P = P(V) = {pij^fj k=i 
with condition (2). Usually [1], [2], [4], [6-11] the matrix P is known. In [3], [5] a constructive 
description of P is given. 

Now we recall this construction. 

Let G = (A, L) be a finite graph without loops and multiple edges, where A is the set of 
vertexes and L is the set of edges of the graph. 

Furthermore, let $ be a finite set, called the set of alleles (in problems of statistical mechanics, 
is called the range of spin). The function a : A — > $ is called a cell (in mechanics it is called 
configuration). Denote by f2 the set of all cells, this set corresponds to E in (1). Let S(A, $) be 
the set of all probability measures defined on the finite set S7. 

Let {Aj,i = 1, ...,m} be the set of maximal connected subgraphs (components) of the graph 
G. For a configuration a G CI denote by cr(M) its "projection" (or "restruction") to M C A : 
a(M) = {a(x)} x< zM ■ Fix two cells a±,a2 G f2, and put 

U(G, <7i, 02) = {cr G SI : cr(Aj) = (Ji(Aj) or <r(Aj) = 02 (Aj) for all i = l,...,m}. 

Now let [i G S(A, <£>) be a probability measure defined on f2 such that //(c) > for any cell 
a e ft; i.e /t is a Gibbs measure with some potential [12], [14]. The heredity coefficients p ai a 2 ,o 
are defined as 



Obviously, p ai a 2 ,<r > 0, Pa lf r 2 ,a = Pa 2 a u <T and Ect6QPctkt 2 ,ct = 1 for all (Ji, cr 2 G fl 

The QSO V = V^ acting on the simplex S^A, $) and determined by coefficients (5) is defined 
as follows: for an arbitrary measure A G S(A, $>), the measure V(X) = A' G S(A, $) is defined 
by the equality 





A '( a ) = Po-io-2,(rA(cri)A(cr2) 



(6) 



for any cell a G 0. 

The QSO construction is also closely related to the graph structure on the set A. 



Theorem 1. [3] The QSO (5) is Volterra if and only if the graph G is connected. 
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3 Non-Volterra QSO 



In this section we describe a condition on measure fi under which the QSO generated by /i 
can be studied using the theory of Volterra QSO. 

Let G = (A, L) be a finite graph and {Aj,z = l,...,m} the set of all maximal connected 
subgraphs of G. Denote by fij = <E> Ai the set of all configurations defined on Aj, i = 1, m. Let 
fii be a probability measure defined on such that //j(cr) > for any <r E fij, i = 1, ...,m. 

Consider probability measure /x on Q = x . . . x Q, m defined as 

m 

i=i 

where a = (o"i, cr m ), with U{ € = 1, m. 

By Theorem 1 if m = 1 then QSO constructed on G is Volterra QSO. 

Theorem 2. T/ie QSO constructed by (5) with measure (7) is reducible to m separate 
Volterra QSOs. 

Proof. Take if = (ipi, <p m ), ip = (ipi, ip m ) G ^- By construction 

Q(G, = jcr = (<Ji,..,<7 m ) G !] : (T; G Y>j}, i = 1, ...,mj 

and 

[ otherwise, 
where we used the following equality 

m m 

H(Sl(G,<p,il>)) = nPi(<7i) = n(jH(<Pi)+tH(il>i))- 

<ri,...,<r m : j =1 j =1 

<T i e{(^ i ,i/' i },i=l,...,m 

Thus QSO (6) generated by measure (7) can be written as 

vw-aw.,,..)- E . fitnff ^ (9> 

Denote 

^i,u; = E = E A(ai,...,CTi_i,W,C7i + i,...,Cr m ). (10) 

CTgQ: o- 1 ,...,CT i _ 1 ,<Tj +1 ,...,<T m 

From (9) we have 



trGfi: trGfi: 



Vl,---,Vi-l,¥i+l,---,>Pm 



n/ij((Tj)l( (Jg^ •,})., .... , , 
— ^ f-^—Xilf!, ...,(pi-i,U,<Pi +1 ,...,(p m )X{lpl,...,lp m ) + 



j 
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y- x 

<Pl,---,<Pm 
*l,---,V'j-l,V'i+l,---,'/'m 

? — '- f-^—X(ipi,...,(f m )X(lp 1 , ^-1,^,^+1, -^m) 

2 y ^ x 

iri,...,(Ti_i,<Tj + i,...,<T m 3=1 rjyr-jj r>j\-rjj 

Note that 

- , fe.- „.. U N(<Pi)+l*M) 

Thus from (11) we have 

RHS of (11) = 

^ A(<^i, ...,<^ i _i,a;,^i + i,...,^ m )A(V'i,...,V'i-i»^)V'i+i) •••,V ; m) + 

(Pli — iVi-l.Vi+l. — .Vm 

2 H 7 VT~^ 7 / a H A(v?i,...,v?j_i,u;,^ + i,...,^ m )A(^i,...,V'm) = 

Vl. — .Vi— liVi+l. — .Vm 



(11) 



_ 2//i(cj) 

ipecii\ui 

Thus operator (9) can be rewritten as 



X - = K + ' (12) 

where Xj i£J is defined by (10), a; € S~2i , i = l,...,m. 

Note that XLeQj = 1 f° r an y * = 1' •••> m - Using this equality from (12) we obtain 

= Xt >- HIh+mS^ • (13) 



V>Gfl 



Comparing this with (4) one can see that for each fixed i (i = l,...,m) the operator (13) is 
Volterra operator V"W : #1^.1-1 _> gl^l-i. The theorem is proved. 

Remark. The set which we used for the set of spin values on G is the same for any 
configuration on fij, z = 1, ...,m. Note that Theorem 2 can be easily extended for more general 
case: when each subgraph Aj has its own <I>j. 
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4 The behavior of the trajectories 



In this section using Volterra QSOs (13) we shall describe the behavior of trajectories of non 
Volterra QSO (9). 



Denote = ^-u^+^-NA ■ It is easy to see that for each fixed i G {1, ...,m} the coefficients 



(i) 

°V i/> satisfy the following properties 

«^ = -«S,l> I°^I<1- (14) 

Thus QSO (13) has the same form with (4). 

If for any i G {1, ...,m} the asymptotical behavior of trajectories of QSO (i.e. (13)) is 
known, say xf^ — > X* u , I — ► oo, then asymptotical behavior of V (i.e. (9)), say A^(cr) — > A* (a), 
Z — > oo, can be found from system of linear equations 

£ A»=X* W , iuen i ,i = l,...,m. (15) 

Now we recall known results ([6], [7]) about trajectories of (4) (i.e. (13) for fixed i). 
Put 

n 

dS 71 ' 1 = {x = (x u ...,x n ) G : Y[ Xj = 0}, int^- 1 = S n ~ 1 \dS n - 1 . 

3=1 

Let x 1 - -* G S 1 " -1 be the initial point. Denote by ^(x°) the set of limit points of the trajectory 
{x^} of the operator (4). 

Theorem 3. [6]. 1) If x^ G intS" 1 " 1 is not a fixed point (i.e. Vx<® / i (0) j, then 

c as™- 1 . 

2) T/ie se£ z/(a;( )) either consists of a single point or is infinite. 

3) If QSO (4) has an isolated fixed point x* G intS" 1-1 , then for any initial point x^ / x* , 
the trajectory {x^} does not converge. 

Corollary 4. The set of limit points v(\ ) of the QSO (9) has properties l)-3) mentioned 
in Theorem 3. 

Proof. 1) By Theorem 3 there is at least one Co G fij with X*^ = 0. Since A(cr) > from 
(15) we get that A* (a) = for all a such that crj = a). This completes the proof of the property 
1). Properties 2), 3) also follow from (15). 

Now we shall recall some notions from theory of tournaments (see [6]) for operator (4). 
Assume ^ for k ^ i. Along with (4) consider the complete graph G n with n vertices. 
Specify a direction on the edges of G n as follows: the edge joining vertices k and % is directed 
from the kth to the ith vertex if < 0, and has the opposite direction if a^i > 0. The directed 
graph thus obtained is called a tournament and denoted by T n . A tournament is said to be strong 
if it is possible to go from any vertex to any other vertex with directions taken into account. 
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Further, after a suitable renumbering of the vertices of T n we can assume that the subtour- 
nament T r contains the first r vertices of T n as its vertices. Obviously, r < n, and r = n if and 
only if T n is a strong tournament. 

Theorem 5. [6]. Suppose that T n is not strong, and let G intS 1 ™ -1 . If j > r, then x^p — > 0, 
ai £/ie rate o/ a geometric progression as / — > oo. 
As corollary of Theorem 5 we have 

Corollary 6. Suppose that there is a io £ {1, ...,m} such that the corresponding tournament 
T^°) of the operator y( l <>) is not strong, and assume X° £ intS 1 " - " 1 . T/ien £/iere is a subset 
£li C f2j smc/i £aai A® (a) — > 0, /or <j£S] wit/i o"j G f2j , af £/ie rate o/ a geometric progression 
as I — > oo. 

Example. Consider graph G = (A, L) with A = {1,2} and L = 0. Take $ = {A, a}. Then 
non-Volterra QSO (9) has the form 

x'i = x\ + 2(i\X\X2 + 2a\X\x^ + 2a\P\X\x^ + 2ai0ix 2 x 3 
x 2 = x 2 + 20 2 xix 2 + 2ai/3 2 x 2 x 3 + 2a\X2Xi± + 2ai/3 2 xiX4 
x 3 = + 2a 2 xix 3 + 2a 2 /3ix 2 x 3 + 20ix 3 x 4 + 2a 2 /3ixix 4 
x' 4 = x\ + 2a 2 2 xix 4 + 2a 2 x 2 x 4 + 20 2 x 3 x 4 + 2a 2 2 x 2 x 3 

where //i = (ai, a 2 ), ay > 0,j = l,2,ai + a 2 = l;/x 2 = (0i,0 2 ),0j > 0,3 = 1,2, Pi + P2 = 1- 

Putting xi + x 2 = -X"i ; i, x 3 + x 4 = Xi >2 and xi + x 3 = A 2j i, x 2 + x 4 = A 2]2 we get the 
Volterra operators like (13): 



X[- L = X hl (l + (2a 1 -l)X 1 ^ 
X' 12 = X h2 (l + {2a 2 -l)X 1 ^ 



(17) 



and 

x 2d = x 2 , 1 (i + m-i)X2,2) 

(18) 



X 2)1 = X 2;1 (l + (2/3 1 -l)X 2;2 ) 
X 2j2 = X 2 , 2 (l + (20 2 -l)X 2;1 ) 



For this example from Corollary 4 and 6 we obtain 

Corollary 7. 1. Any trajectory of non-Volterra QSO (16) has following limit 

{ (1,0,0,0), if 2ai > l,2/3i > 1, 
(0,1,0,0), if 2ai > l,2/9i < 1, 
1^00 ) (0,0,1,0), if 2ai < l,2/3i > 1, 

(0,0,0,1), if 2ai < l,20i < 1. 
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2. If 2Pi = 1 then S\ = {x : X3 = X4 = 0} and S2 = {x : x\ = X2 = 0} are the sets of fixed 



3. // 2ai = 1 then S3 = {x : X2 = x^ = 0} and S4 = {x : x\ = X3 = 0} are the sets of fixed 



4. If 2a\ = 2b\ = 1 then S5 = {x : X2 = x^x\ = X3} and Sq = {x : x\ = X2,X3 = X4} are 
i/ie sete of fixed points for (16). 
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